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ABSTRACT 



The hydrodynamic stability of an underwater station was studied 
using wind tunnel models. Rigid model force and moment coefficients 
were obtained for various orientations of the station in three 
dimensions. The unsteady motion of the station in cross currents 
was investigated with a dynamically similar model. Color 8 mm 
pictures and amplitude- time traces of the unsteady motion were 
obtained. 
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I . INTRODUCTION 



The purpose of this study was to investigate the hydrodynamic 
stability of a manned underwater station when the station is fully 
submerged and subject to cross currents. The concept of a manned 
underwater station, or MUS, has been investigated by General Dynamics 
Corporation, Electric Boat Division, under contract for the U. S, 

Naval Civil Engineering Laboratory, Port Hueneme, California (Ref. 1). 
The present study, which is in support of the earlier MUS work, was 
sponsored by NCEL under work order WR-8-0128 and was monitored by 
Dr. Cheng Lung Liu of NCEL Ocean Engineering Department - 

The object of the MUS project as stated in Ref. 1 is to develop 
a deep ocean submergence habitat which will accommodate a five man 
crew for an operational period of thirty days. The station is to 
be capable of submergence to depths of 6,000 feet. 

As envisioned by General Dynamics, the station would consist of 
two circular cylinders of unequal diameters and heights. The smaller 
diameter cylinder is capped with hemispherical ends and will contain 
a nuclear reactor. The larger diameter cylinder is capped with 
spherical ends which are fitted to the cylinder by conical, tangential 
shells. (Fig. 1) The habitat cylinder will have five decks and is 
the working space for the crew. 

Since these cylinders are designed only for a vertical orientation, 
a large ballast tank completely encircles both cylinders. The design 
of this tank is based in part on weight and balance considerations. 
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The station is not equipped with a propulsion device although 
it does have ballast and buoyancy control equipment . The concept 
of the station is somewhat unique and the MUS may best be described 
as an elevator. In concept, the station is to be towed on the surface 
to a pre-selected site by a support vessel. When the operations 
site is reached, the station will be manned and an emplacement 
procedure will commence. The MUS will utilize an anchor and cable 
arrangement to pull itself to the bottom in bootstrap fashion. It 
is equipped with a winch (located on the main ballast tank) and 
attendant systems for control. 

When submerged, the MUS will maintain a slight positive 
buoyancy at all times. The anchor will first be lowered to the 
ocean bottom while the MUS holds on the surface. Once the anchor 
is in place, the buoyancy of the MUS is reduced and winchdown begins. 
As the winch retrieves the cable, the MUS will haul itself to the 
bottom. The recovery cycle, i.e., the ascent to the surface, is 
just the reverse of this procedure. 

It is primarily with the emplacement and recovery operations that 
this study is concerned. When the MUS is in the mid-water column, 
it may well be subjected to a cross flow caused by deep ocean 
currents. Evidence exists to expect such currents throughout the 
water column, having a magnitude of 0,1 to 0.5 knots. 

If such currents were to exist, it would be reasonable to expect 
the MUS to be forced into motion until equilibrium is attained. 

Such motion could be translational or rotational and perhaps 
oscillatory in nature. Since the MUS is restrained only by a single 
cable (to the anchor), the only other stabilizing elements are: 
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1. The buoyancy couple caused by the displacement of the meta 
center relative to the center of gravity. 

2. Hydrodynamic stabilization, i.e., the geometry of the station, 
including such extremities as the legs and pads. 

3. The inertial resistance to motion (moment of inertia) 

Another aspect of the MUS recovery cycle which necessitated a 

stability investigation was that of emergency ascent. It was proposed 
that, in the event of an emergency such as cable jamming or bottom 
entrapment of the support pads, a free ascent could be accomplished 
by cutting the cable (with explosive cutters) and/or jettisoning the 
support pads and legs. Since these are in fact stabilizing devices, 
such an ascent could only be feasible if the resulting configuration 
of the MUS provided sufficient stabilization. 

The objectives of the experiments reported herein were twofold, 
namely: 

a. To establish by means of rigid model wind-tunnel tests the 
pertinent dimensionless force and moment coefficients for 
the MUS at various orientations, and 

b. To determine the dynamic behavior using a f, partial mode" 
dynamic model, and in particular to note the existence of 
any self-induced and/or vortex-induced oscillations. 

In this investigation, the MUS was represented by an 0.0209 
scale model. Reynolds numbers were matched between model and prototype 
so that aerodynamic measurements on the scale model were applicable 
to the hydrodynamic traits of the prototype. 

One should note that the MUS configuration is still in the early 
design stages, hence the results obtained in this study should have 
as a goal the possible application to other configurations. However, 
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accurate forecasting tends to be uncertain because the bluff-body 
features of these configurations do not interact in a linear manner. 
The testing accomplished in this study was relatively easy, and other 
configurations could be studied using a similar approach. 
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II DISCUSSION 



A. DESCRIPTION OF MODELS AND EXPERIMENTAL PROCEDURE 

The investigations described herein upon the 0.0209 scale MUS 
models were performed at the Naval Postgraduate School (NPS) low- 
speed wind tunnel facilities. The force and moment measurements 
were conducted on the "rigid” model in the Aerolab 32 in. x 45 in. 
wind tunnel while the scaled dynamic model was tested in the West 
Coast Research 48 in. x 60 in. wind tunnel. 

The "rigid" model was linearly scaled and was fabricated from an 
overlay of Honduras mahoganyy upon an aluminum spar or reinforcing 
plate (Fig. 6). The purpose of the spar plate was to make possible 
the attachment of the support arms and tail sting at prelocated or 
orientations of the model. It was possible to position the model 
in a multitude of ways upon the conventional three strut support 
system so that representative yaw and pitch attitudes could be 

+ ° 4. o 

obtained for an angle range of -180 and -90 respectively. Figure 
5 presents a sketch of the model including pertinent dimensional 
data relative to the force model, while Figures 11 to 44 show the 
model installed in the wind tunnel at various orientations. The 
model was tested at dynamic pressures from 3.85 to 57.8 psf (in air) 
which correspond to prototype velocities in sea water from 0.08 to 
0.30 knots (refer to Figure 52 for conversion from model to prototype 
operating conditions) . 

The forces and moments were measured upon an Aerolab three- 
component beam balance (Fig. 10) which had the capability of measuring 
lift, drag, and pitching moment relative to the wind axes. 
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The balance limits were: 



Component 


Limits 


Sensitivity 


Lift force 


- 150 lbs 


- 0,010 lbs 


Drag force 


- 50 lbs 


- 0.005 lbs 


Pitch moment 


- 1000 in-lbs 


- 0.050 in-lbs 



Care should be exercised in interpreting the meaning of lift and 
drag forces and pitching moments at some of the extreme model positions, 
such as when the model has been rolled 90 degrees on the balance, 
since in this attitude pitch angle changes actually corresponded to 
model yaw angle changes (Figs. 43 and 44). In this extreme case, 
lift force and pitching moment correspond to side force and yawing 
moment respectively. It is unfortunate that the balance system was 
not readily usable as a six-component type since then it would have 
been possible to convert the wind axes results to stability or body 
axes as is frequently done in aircraft analysis. 

Data corrections were made for weight tares, while aerodynamic 
tares to account for the influence of the support system were 
approximately estimated by measuring the balance force and moments 
without a model present. Since the model was a relatively bluff body, 
as compared to a streamlined aircraft configuration, these 
approximations were considered as being reasonable. 

The dynamic model tests may be described as "partial mode" tests 
since only selected degrees of freedom were investigated. The modes 
investigated included lateral and axial translation and yaw rotation. 

The tests were performed upon two 0.0209 scale models, one made from 
balsa and the other from sugar pine (Figs. 7, 8, and 9). 
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The balsa model had as its purpose an initial observation o£ the 
models dynamic behavior and an experimental confirmation of the 
yawing moment static trim points. The full scale weight and weight 
distribution table for the MUS from Ref. 1 was used as a basis for 
design of the pine model. By using a component approach, the full 
scale weights and weight distribution were scaled to obtain the desired 
mass moment of inertia for the pine model. Removable lead inserts 
were utilized to allow a variation of mass and mass moment of inertia. 
The balance of the model about the desired center of gravity location 
was checked experimentally, while the mass moment of inertia was 
determined by utilizing the principle of the torsional pendulum. 

An object having a known moment of inertia, namely a solid steel 
cylinder of given weight and geometry, was suspended on an 0.037 
inch diameter steel piano wire of known length. The oscillation 
frequency was then measured which in turn allowed the determination 
of the wire's shear modulus of rigidity. From the physical properties 
of the wire (diameter, length, and shear modulus), and the measured 
free vibration oscillation frequencies, it was possible to determine 
the mass moment of inertia of the model. Table II lists the physical 
properties of the dynamic models tested. 

A narrow diameter tube was inserted along the vertical axis 
(through the assumed prototype center of gravity) in order to shield 
the dynamic model from the support wire. The model was suspended on 
an 0.037 inch diameter steel piano wire which was fixed to supports 
external to the tunnel. (Figs. 56 through 59) The wire was preloaded 
to 50 lbs. tension by means of turnbuckles, and the tension was 
frequently checked by using a standard rigging tensiometer. Initially, 
only the taut wire support was used, but later tests were conducted 
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with the translational degrees of freedom removed by means of local 
tie-down constraints. In addition, the last sequence of dynamic 
model tests were provided with an electrical angle sensing potentiometer 
in order to allow time-history traces to be made of the angular motion, 
(Figs, 60 and 61) 

Motion picture photography of the dynamic model tests was accom- 
plished by using a mirror arrangement. An adjustable mirror was 
mounted to the tunnel wall above the model so that overhead photography 
was possible. The later tests involved simultaneous photography 
and electrical data collecting, (Table VI) 

B. SCALING CONSIDERATIONS 

The scale of the models (0,0209) was established on the basis of 
the largest model compatible with the wind tunnel test section size 
and velocity range in combination with a Reynolds number match 
between model and prototype. 

For the prototype design with an ocean current of 0 P 3 knots in 
salt water, an estimate of full scale Reynolds number would be 



The above Reynolds number is based upon a mean diameter of 14,93 
feet. An 0.0209 scale model ( X = 0,0209 = scale ratio) would satisfy 
the prototype Reynolds number condition above at a tunnel velocity of 




N Re = 4. 4-5“ x /O 



(prototype in salt water) 



V 
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which corresponds to a tunnel dynamic pressure of 60.3 psf. The ratio 
of model to prototype velocities is 443 under the assumption of 
equivalent Reynolds numbers. 

The mass moment of inertia scaling of the dynamic model was 
established by the free vibration equation in the yaw mode. As 
an approximation, one may consider the slope of the C N vs V curve 
to be linear in the immediate neighborhood of a zero moment static 
trim position. 

The total yawing moment may be approximated by the sum of the 
linearized static term and a term due to the dimensionless angle 

rate (%rV 

c N -- L c » r 1 v * Tv t (1) 

Multiplying ^ by to obtain physical units, then the free 

vibration equation in yaw is: 



■±^11 = ° 



dt 



(2) 



where : 



= Mass moment of inertia, taken about the £ axis. 
Dividing through by — g^- an< ^ introducing the dimensionless time, ft 



where: “ft — b 



(-SL\ 



then 



/2y\ 4^ _ C N 41 _ c Y = o 

''d 1=* y it- f ' 



(3a) 



The enforcement of Reynolds number equivalence insures that the 
dimensionless yawing moment coefficient i s identical for 

both model and prototype. In order to achieve full dynamic similarity, 
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the nondimens ional , natural frequency of the model must be identical 
with that of the prototype. This argument then defines the scaling 
law for the mass moment of inertia as shown below. 



at 1 






= o 



(3b) 



This may be compared to the familiar vibration equation 



Kit 2 ' 



^2? CO, 



41 

dt- 



+- 60 h 1=0 



(4) 



In order for the equation of motion to describe concurrently the 
behavior of both model and prototype, each term must equate in 
model to prototype conversion, therefore 




vnod ie\ 




Pro+o^pe 



or 









Wv 






gince 



f' J p 



then 



=( C V> 

= TO 

^ X 2 - 

Ar 

^ = A 



a_ 

d -1 



t . A 3 



it 



m 







(5) 



if the inertial scaling of above were used, then the non-dimensional 
frequency would be alike between model and prototype, i.e., 
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Similarly, a check of the damping terms: 



>/\ ? 

L-i* 







will yield the same result by the same argument. 



Unfortunately, such scaling was not physically realizable as it 
would demand a prohibitively light construction of the model. 

Returning to the non-dimensional natural frequency, which may 
be identified in eqn. (3b); 



- 2 






1 r 

8 X* 



(4a) 



_ ( ^>L \ )\ ( \ 

CCd^f K J V- ) 



if the mass moments of inertia were set as: 







then the dimensionless frequencies won’t correspond, but will scale 

(fin _ P ^ _ i 2.3 x 10 3 

(®h* D p f f 

Considering that dimensionless frequency is by definition: 

i_ tl , 

' 2v 



one readily obtains that: 







( 6 ) 
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After substitution of the scaling parameters (cf. Table I), the true 
frequency scaling becomes: 



Cw^p = o. 00 134 

For example, if ~ 0,1 0 CpS } then 

= 0.001-54(0.70 ^(60 ^ ;n ) = 0,o56> 2 C ^j~ S 



and the prototype period would be 18.1 min/cycle. 

Using a quasi-steady argument, enforced by the case where the 
damping ratio ^ is zero, one may use the above scaling law 
to convert model frequencies to prototype frequencies for the yaw mode. 

In a similar fashion to the above analysis concerning the yaw 
mode, one may consider the translational mode free vibration equation. 
In the lateral translational mode, a similar analysis to the yaw 
mode scaling may be made starting from the free vibration equation, 

^ £ ' C F^ (V) + $ti ft = o 



(7) 



The damping derivative, Cp . , may be approximated by the quasi- 

% 

steady assumption in a manner similar to that used by Parkinson, 

(Ref. 2) when he considered the limit-cycle solution for the galloping 
transmission line problem. Using the sign convention of Figure 2, 
one obtains that: 



C F • ~ + v 



y 



since 



+ k 



Substitution of this approximation into equation (7) yields that: 



u - it 



% 



mV 7^ 



C = + •fe % = ° 



( 8 ) 
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It is convenient to convert equation (8) to dimensionless time, £ , 

'V Ui j 

and dimensionless lateral translation, Q , where ^ — /d 

Direct substitution plus a slight amount of algebraic simplification 

yields : 



fOL - C F . d£ + 4 £ 






4^ 'r d£- 



m 4 v 



^ ^ = o 



(9) 



Equation (9) may be compared with the second order differential 
equation; equation (4), in order to identify that in the lateral 
translational mode, the non-dimensional frequency is: 




where _ ci 




For the model installation, the spring constant, k, is proportional 
to the tension in taut wire support as may be seen in the rear 
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An approximation for the tethered prototype vehicle may be made on 
the assumption that the cable tension is equivalent to the buoyancy 
force. For this case, one would obtain 



where the tether cable length, JL p , is an equivalent value which 
may be estimated. 

The initial consideration will be to find the true frequency value 
of the prototype on the assumption of equivalent non-dimensional 
frequencies, followed by a correction for differences in cable tension. 
For the case of / — v / -r 

(iO p = C^n)* 



it follows that: 
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Next, we shall investigate the implications of cable tension. Since 
, it follows that: 



/jfep xt = i ms. 



m Y ( 



G> JL 



(Vjy 2 ev( 






or 



= 



46oo 



*145-0 J 



-\ 



67, 8 






If cable lengbbs are approximately given by Xp = 600 ft., 
J? - 6 ft. and model cable tension is 50 lbs., then 



1 P 

su 



2 . 1)00 



hn 



X 






or 



r p = 2,ie ii 
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Actually it is envisioned that the prototype buoyancy force will be 
on the order of 100 lbs., hence the corrected frequency would be 

Corf ' 

~ 3/70 

Therefore if the lateral translational frequency of the model were 
about 1 cps , then the representative prototype period would be 
about 53 minutes per cycle. 

In retrospect, the dynamic scaling problem would require 
simultaneous scaling of mass, mass moment of inertia about the 
yaw axis, and cable tension. Since this procedure was not initially 
followed, estimates were made to account for the relative effects 
of the mismatches. 
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Ill, RESULTS 



A, RIGID MODEL DATA 

Various rigid model plots are included (Figs. 45 through 55). 

Test runs 1 through 37 indicated that the hydrodynamic coefficients 
are essentially constant throughout the Reynolds number range under 
consideration. A large number of Coef f icient-Reynolds number plots 
were obtained (see run schedule. Table VI), and representative 
plots are shown in Figures 46 through 48 inclusive. The largest 
values of the hydrodynamic coefficients that were measured, with the 
corresponding full scale forces and moments, are plotted for 
comparative purposes in Figure 52. These forces, as compared to 
the buoyancy force of 100 lbs., are considered to be significant. 

A drag coefficient contour plot is included in Figure 45. The 
drag coefficient values along the ordinate and absisca of this plot 
were measured in five-degree increments (See also Figs. 49 and 54), 
while other values were measured in combinations of thirty-degree 
increments . 

The plot of yawing moment coefficient versus yaw angle (Fig. 53) 
shows that the MUS has two stable trim points in yaw, at zero degrees 
and 180 degrees, while two unstable trim points exist at 90 and 270 
degrees of yaw. It should be emphasized that these are zero yawing 
moment positions, while the slope indicates the degree of stability. 
Also note that stable trim did not occur at the orientation for 
minimum drag. 

Direct measurement of the yawing moment vs yaw angle has 
application to defining the aerodynamic restoring moment term in the 
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yaw-degree-of-f reedom dynamic equation. However, means do not exist 
for direct measurement of the damping term due to yaw rate.. The 
damping term must be deduced indirectly by observing the time history 
of decaying yaw oscillation. 

The plot of side force versus yaw angle at zero pitch angle 
is shown on Figure 55. An approximate analytic expression for the 
variation of side force coefficient is: 



However, the actual plot shows many irregularities in the curve 
which may be related to the complex flow interactions between the 
various bluff-body components that comprise the MUS. Especially 



It should also be noted that the above expression for side force 
coefficient will result in the damping term of the lateral trans- 
lational degree of freedom having a negative value. The consequence 
of this feature was noted in the dynamic model as a limit cycle 
type of motion. 

B, DYNAMIC MODEL RESULTS 

A color 8 mm motion picture film was prepared as a supplement 
to this report* Part one of the film shows the experiments with 
the single wire arrangement, and careful observation will show the 
lateral motion beating phenomenon. In this configuration the model 
had degrees of freedom both in yaw and lateral translation. The 
combination of an apparently lightly damped yaw mode plus a 
negatively damped lateral translation mode produced a complex 
response with an apparent indication of beating due to the closeness 



c 



F 
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noteworthy is the rapid unstable trim points, 
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of the two frequencies. The lateral excitation caused the self- 
induced nature of the oscillation in yaw, but this did not appear 
when the lateral constraint was imposed. 

Part two of the film shows the instrumented and constrained 
model experiments. No lateral motion occurred and the damping 
ratio was small but finite until a dynamic pressure of five or six 
psf was reached. Strip recorder charts are included. (Fig. 67) 
Dynamic model data are summarized below in Figs. 62 through 66. 

The predicted frequencies for the model, shown as estimates on 
Figs. 62 to 64, are based upon the recognition of frequency from 
equations (2) and (4a) as: 




Although the predicted frequencies are about 10 to 15 percent higher 

i 

than measured values, the trends with dynamic pressure are proper 
and add confirmation to the concept of an equivalent hydrodynamic 
torsional spring term. It is quite possible that the actual observed 
frequencies reflect rate dependence of the hydrodynamic torsional 
spring term on the bluff body. 
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IV, ANALYSIS 



A. STATIC STABILITY 

The MUS stability in the pitch mode may be interpreted from Figure 
50; i.e., the pitching moment coefficient versus pitch angle plot. 

This plot was obtained only for zero degrees of yaw. At minus ninety 
degrees pitch angle, i.e., the ascending mode, the MUS will trim 
but be statically unstable. While the pitching moment itself is 
zero at o( = -90 , any increase in pitch angle results in a positive 
pitching moment which tends to increase the pitch angle further. 
Similarly, a decrease in pitch angle results in a negative pitching 
moment. At plus ninety degrees pitch angle i.e. the descending mode, 
the MUS is again statically unstable. In the neighborhood of zero 
degrees pitch angle, the MUS has a positive pitching moment coef- 
ficient and hence is not trimmed. Therfore, in the absence of the 
cable constraint, the MUS is statically unstable in the pitch mode. 

The balsa wood dynamic model was suspended in a manner so as to allow 
a pitch degree of freedom for one set of exploratory investigations. 

Its dynamic behavior tended to confirm the preceding conjectures since 
the stable trim points appeared to be weak and the model was quite 
prone to rotate in pitch. 

In the vertical translation mode the MUS is statically unstable. 
Referring to Figure 51^ =0, ^ = 0 , , in a direct cross current, 

the MUS has a negative lift, which will tend to cause the MUS to 
descend. This in turn would increase the effective pitch angle, 
and result in a more negative lift. 
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Figure 53 shows that the station is statically stable in the yaw 
mode (at zero pitch) at yaw angles of zero and 180 degrees and 
statically unstable in yaw at 90 and 270 degrees of yaw. 

Figure 55 shows that the station is statically unstable in lateral 
translation at yaw angles of zero and 180 degrees. If a lateral 
translational velocity occurs, an effective yaw angle is introduced 
in the nature of sideslip. When this occurs, the side force is 
such that the lateral translation tends to be further increased 
(forced) . 

B. DYNAMIC STABILITY 

In the absence of constraints, the experiments indicated that the 
MUS would be dynamically unstable in all translation modes and in 
the pitch mode. 

The yaw mode appeared to be dynamically stable, as the oscillation 
amplitude either decreased in time or remained constant. Negative 
damping ratios were not observed for the range of tunnel speeds used, 
but it is possible that higher dynamic pressures could have induced 
negative damping. 

The scaling relations for the yaw mode and the lateral translation 
mode indicate that the natural frequencies of these modes are very 
widely separated at full scale. This would indicate that the coupling 
tendency observed in the model experiments would not occur for the 
full scale MUS. Thes most important conclusion from scaling 
considerations is that an oscillation on the model in a wind tunnel 
on the order of a second per cycle period would scale to the prototype 
vehicle in water with a period on the order of 20 to 50 minutes per 
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cycle. Therefore it is quite possible that unusual oscillational 
behaviors in the tunnel would not be perceptible to the crew enclosed 
in the MUS because of the low velocities involved. 

C. COMMENTS ON MUS CONSTRAINTS 

The experimental results indicated that the proposed free ascent 

mode would be practicable. Considering an assumed stabilized condition 

where the MUS would be floating freely in any existing current, with 

cable and pads jettisoned, and in the absence of an unequilibrated 

buoyancy force, no hydrodynamic force would be generated. If the 

positive buoyancy force of approximately 100 lbs. were considered, 

the MUS would be forced into vertical translation and hydrodynamic 

forces would be generated. In this orientation, the wind tunnel 

axis position of minus ninety degrees of pitch angle would be 

applicable (Fig. 2). Assuming that a constant side current existed 

throughout the water column, the MUS would begin to rise, sensing 

. o 

only the hydrodynamic forces corresponding with o( = -90 . 

However, in this orientation, both a drag force and a lift force could 
be generated. If one assumes that the MUS could achieve a vertical 
velocity of 0.2 knots, from the wind tunnel results, the resulting 
drag force would be 100 lbs., and the lift force would be zero 
(Figs. 49, 51, 52), Actually, the resulting drag forces would limit 
vertical ascent speed to 0.2 knots, if the 100 lb. buoyancy force were 
maintained. This speed limitation may be verified by cross referencing 
Figs. 49 and 52. The drag coefficient of 0.7 will result in a drag 
force of 100 lbs. at 0.2 knots, therefore the drag force will 
equilibrate the net buoyancy force at 0.2 knots. A check of the 
pitching moment vs. pitch angle plot (Fig. 50) for this orientation 
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shows that the MUS would be at an unstable trim point, that is the 

pitching moment would be zero but the variation of moment with flow 

angle would not produce a restoring moment. As long as the flow 

j o 

vector remains constant at 0\ = -90 , the MUS will remain upright. 

If the flow vector is disturbed from 0( = -90° or more likely, if 
the MUS begins ascending from an off-vertical orientation, the 
resulting pitching moment will be in a direction tending to over- 
turn the MUS. This pitching moment would be on the order of 300 

o 

ft-lbs. maximum (for a 0.2 knot vertical speed and a 20 off-vertical 
orientation). Considering the design (buoyancy couple) stability, 
this slight pitching moment should present no difficulties (ref. 1). 

That is, the location of the center of buoyancy above the center 
of gravity results in a stabilizing moment. 

Considering the case where the MUS is tethered by the cable in 
the mid-water column, if side currents exist with a speed of 0.3 
knots, a pitching moment would be generated on the order of 1300 
ft-lbs. maximum. This would place a side load on the cable attach- 
ment point, as the cable would equilibrate the pitching moment. In 
addition, a yawing moment would exist, having the same order of 
magnitude. The yawing moment would apply a torsional load to the 
cable or to the cable attachment point, depending on the actual 
design, however, the vehicle would tend to align itself to the 

/ o o 

stable trim point corresponding to lu = 0 or 180. 

The oscillatory behavior of the model appears to be of no 
practical consequence to the full scale MUS and probably would not be 
noticable by the crew. Considering the proposed emplacement procedure, 
since the principal static equilibrating force that exists to counteract 
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drag forces (in side currents) is the 100 lbs. buoyancy force, the 
wind tunnel measurements indicate that a forced descent would result 
with current speeds as small as 0.2 knots. At this speed, the resulting 
drag force would be equal to the buoyancy force, and therefore, static 
equilibration would not be possible. 

The restoring force is equal to: 

F = T sin 

X 

Where T - buoyancy force = cable tension 

and ^ = deflection angle of the cable 

from the vertical 

if T = 100 lbs. 

for a current speed of 0.15 knots and = 0.7 

D = 60 lbs. (Fig. 52) 

Sin ^ = 0.6 

= 37° 

or in other terms, a horizontal translation of 3,600 ft. and a 1,200 
ft. descent would occur before static equilibration if 6,000 ft. of 
cable separated the MUS and its anchor. For this reason, it is 
recommended that a higher buoyancy force be considered. 

The foregoing discussion on forced descent is admittedly overly 
simplified and does not take into consideration the special effect 
of catenary behavior of the cable itself. This topic would merit 
individual study and analysis before a complete discussion of 
forced descent could be made. It is also recognized that increasing 
the buoyancy force will create a serious problem involving free 
ascent terminal velocity. In this respect, a compromise would be 
involved between stability considerations • and ascent problems. 
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TABLE I 



VALUES OF SOME PHYSICAL CONSTANTS 



Characteristic area (model). • . . 
Characteristic area (prototype) 

Height of model 

Height of prototype ....... 

Characteristic length (model) * . 

* 

Characteristic length (prototype) 
Scale ratio, X 



X 1 




-u, salt water 
"V, air . . . 



Dimensionless velocity, V 




. . . 0.546 ft 2 

. . . 1248 ft 2 

. . . 1.0 ft 

. . . 47.8 ft 

. . . 3.75 in 

. . . 14.93 ft 

. . . 0.0209 

. . 4.37 x 10 4 

. . 9.0 x 10" 6 

. . 4.0 x 10' 9 

0.17 x 10 4 ft 2 /sec 
1.58 x 10 4 ft 2 /sec 
. . 1.23 x 10~ 3 

. . . 443 



Note: The characteristic length was arbitrarily chosen to be the 

average or mean of the diameters of the five principal 
components which comprise the MUS. 
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TABLE II 
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TABLE IV: LIFT COEFFICIENT S 
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TABLE VI: RUN SCHEDULE 



A. Rigid Model Tests, Model A in Aerolab Wind Tunnel 



Dynamic 

Test Run Pressure, 

Number (psf) 



Pitch 

Angle 



Yaw 

Angle 



Comments 



1 


3.85 to 57.8 
in increments 
of 4.0 


0 


0 


2 


II 


-30 


0 


3 


II 


+30 


0 


4 


II 


0 


30 


5 


It 


-30 


30 


6 


II 


+30 


30 


7 


II 


0 


60 


8 


II 


-30 


60 


9 


II 


+30 


60 


10 


II 


0 


90 


11 


It 


-30 


90 


12 


It 


+30 


90 


13 


II 


0 


120 


14 


II 


-30 


120 


15 


II 


+30 


120 


16 


II 


-30 


150 


17 


II 


0 


150 


18 


II 


+30 


150 


19 


1 1 


0 


180 


20 


II 


-30 


180 



C, c vs N 



fc •> L j'-M 



Re 
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21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 



Dynamic 

Pressure, Pitch 


Yaw 




(psf) Angle 


Angle 


Comments 


3.85 to 57.8 +30 

in increments 


180. 


N RE 


of 4.0 






" 0 


210 


II 


" -30 


210 


II 


" +30 


210 


II 


" -30 


240 


II 


" 0 


240 


II 


" +30 


240 


II 


" 0 


270 


A vs K* 


" -30 


270 


II 


" +30 


270 


II 


" 0 


300 


II 


" -30 


300 


It 


" +30 


300 


II 


" 0 


330 


II 


" -30 


330 


II 


" +30 


330 


II 


46.1 -30 to +30 in 

5 increments 


0 


C 0 A C M VS o( 


3.85 to 57.8 +90 

in increments 


0 


C o ) C l,C m VS h/ R 


of 4.0 






46.1 +60 to +120 in 


0 




5 increments 






3.85 to 57.8 +60 

in increments 

n-P /. r\ 


180 


VS 



44 



41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 



Dynamic 
Pressure , 


Pitch 


Yaw 




(psf) 


Angle 


Angle 


Comments 


3.85 to 57.8 
in increments 


+60 


0 


C o, C L C M vs N 


of 4.0 

tl 


+60 


30 


II 


tl 


+60 


210 


11 


34.62 to 46.1 


+90 


30 


II 


3.85 to 42.3 
in increments 


+60 


60 


II 


of 4.0 


+60 


240 


II 


30.8 to 38.5 


+90 


60 


II 


3.85 to 42.3 


+60 


90 


II 


3.85 to 38.5 


+90 


90 


II 


3.85 to 42.4 
in increments 


+60 


270 


° D c C M VS IV 

^ i L ) M 


of 4.0 

11 


+60 


120 


it 


II 


+60 


300 


it 


II 


+60 


150 


ii 


11 


+60 


330 


ii 


46.1 

in 


-120 to -60 
5 increments 


0 


,c„ vs 


3.85 to 46.1 
in increments 


-60 


180 


C »A 


of 4.0 








It 


-90 


0 


II 


11 


-60 


0 


II 


11 


-60 


210 


11 



45 



60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 



Dynamic 

Pressure 

(psf) 

46.1 

II 

II 

II 

It 

II 

It 

It 

II 

II 

II 

II 



Pitch Yaw 

Angle Angle Comments 

0 345 to 015 in C Q Q yc, llr 

5 increments ) P ) h/ 7 

" 020 to 045 in " 

5 increments 

" 050 to 075 in " 

5 increments 

" 080 to 105 in " 

5 increments 

" 110 to 135 in 

5 increments 

" 140 to 165 in " 

5 increments 

" 170 to 195 in 

5 increments 

" 200 to 225 in 

5 increments 

" 230 to 255 in 

5 increments 

" 260 to 285 in 

5 increments 

" 290 to 315 in " 

5 increments 

" 320 to 340 in 

5 increments 
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B. Dynamic 


Model Tests, West 


Coast Research Wind 


Tunnel 


Test Run 
No. 


Model 

Used 


Dynamic 
Pressure , 


Photographs 
(Motion Picture) 


Time-History 

Trace 


Not Numbered 
(Note 1) 


B 


0.5 to 
l'.Q'pSf 


none 


none 


M-l - M-15 
(Note 2) 


C 


5.0 


yes 


none 


1-9 

(Note 3) 


C 


5.0 


no 


yes 


10, 11 


c 


2.5 


no 


yes 


12 


c 


7.5 


no 


yes 


13, 14 


c 


2.0 


no 


yes 


15, 16 


c 


3.0 


no 


yes 


17, 18 


c 


4.0 


no 


yes 


19, 20 


c 


5.0 


no 


yes 


21, 22 


c 


6.0 


no 


yes 


23 


c 


7.0 


no 


yes 


24, 25 


c 


8.0 


no 


yes 


24A 

(Note 4) 


c 


2.0 


yes 

film test 
sequency 16 


yes 

(note 5) 


25k 


c 


3.0 


yes 

film test 
sequence 17 


yes 


26 


c 


4.0 


yes 

film test 
sequence 18 


yes 


27 


c 


5.0 


yes 


yes 
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Test Run 
No. 


Model 

Used 


Dynamic 
Pressure , 


Photographs 
(Motion Picture) 


Time-History 

Trace 


28 


C 


5.0 


yes 

film test 
sequence 19 


yes 


29 


C 


6.0 


yes 

film test 
sequence 20 


yes 


30 


c 


7.0 


yes 

film test 
sequence 21 


yes 


31 


c 


7.0 


yes 

film test 
sequence 22 


yes 


32 


c 


8.0 


yes 

film-test 
c sequence' 23 


yes 


33 


c 


8.0 


yes 


yes 


L-l , L-2 , 
L-3 


D 


2.0 


no 


yes 


L-4 , L-5 


D 


3.0 


no 


yes 


L-6 , L-7 


D 


4.0 


no 


yes 


L-8 , L-9 


D 


5.0 


no 


yes 


L-10, L-ll 
L-12 


D 


6.0 


no 


yes 


L-13, L-14 


D 


7.0 


no 


yes 


L-15 , L-16 
L-17 , L-18 


D 


8.0 


no 


yes 


LL-1 , LL-2 


E 


2.0 


no 


yes 


LL-3, LL-4 


E 


3.0 


no 


yes 


LL-5 , LL-6 


E 


4.0 


no 


yes 


LL-7 , LL-8 
LL-9 


E 


5.0 


no 


yes 
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Test Run 
No. 


Model 

Used 


Dynamic 
Pressure , 


Photographs 
(Motion Picture) 


Time-History 

Trace 


LL-10, 


LL-11 


E 


6.0 


no 


yes 


LL-12 , 
LL-14 , 


LL-13 

LL-15 


E 


7.0 


no 


yes 


LL-16, 


LL-17 


E 


8.0 


no 


yes 



Notes : 

1. Various runs conducted in both the Aerolab tunnel and the 
West Coast Research tunnel. Model mounted on wire in order 
to observe both the yaw mode and the pitch mode. 

2. Local tie down constraints not used; these runs correspond 
to Part One of the movie film. Film test sequences 1 
through 15; film speed, 16 frames per second. 

3. Local tie down constraints were used in this and subsequent 
runs listed. 

4. Runs 24A through 33 correspond to Part Two of the movie 
film. Film test sequences 16 through 23; film speed, 

16 frames per second. 

5. Test run numbers (24A through 33, L-l through L-18, LL-1 
through LL-17) correspond with notation used on the 
time-history traces . 
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FIG 1 MOCK-UP OF THE MANNED UNDERWATER STATION 





2. FORCES AND MOMENTS ARE RELATIVE TO 
WIND AXIS SYSTEM (AS ABOVE). 



FIG. 2 COORDINATE AXIS SYSTEM 
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z 




FIG. 3 YAW ANGLE DEFINITION 
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FIG. 4 PITCH ANGLE DEFINITION 
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FIG. 5 

EXTERIOR DIMENSIONS 
OF THE MUS MODEL 
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FIG. 6 ALUMINUM SPAR 

RIGID MO 
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FIG. 7 

WEIGHT DISTRIBUTION 
IN THE DYNAMIC MODEL 




WEIGHTS A , B , C , D ,E ,F, 
AND l - 1.0 In. dia 
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COMPONENT 


WEIGHT 


DIA. 


HEIGH T 


WEIGHT A 


0.30 


lbs 


1.0 


0.96 


WEIGHT B 


0.G7 


lbs 


1.0 


207 


WEIGHT C 


0-67 


lbs 


1.0 


2.07 


WEIGHT D 


0.67 


lbs 


1.0 


2.07 


WEIGHT E 


0.67 


lbs 


1.0 


2.07 


WEIGHT F 


0.67 


lbs 


1.0 


2.04 


WEIGHT G 


1.45 


I bs 


2.0 


1.12 


WEIGHT H 


1.00 


l bs 


2.0 


0.7 7 


WEIGHT 1 


0.31 


lbs 


1.0 


0.96 


HOLDING 

TANK 


0.20 


lbs 






REACTOR 

CYLINTER 


0.74 


lbs 






HABITAT 

CYLINDER 


O.PC 


lbs 






MAIN BA. AST 

TANK 


1.25 


lbs 






OBSERVATION 

SPHERE 


0.05 


lbs 







FIG. 8 

WEIGHT DISTRIBUTION 
IN THE DYNAMIC MODEL 
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FIG. 10 THE AEROLAB 3 COMPONENT BEAM BALANCE 
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8 min 
13* max 




Fig* 11 - Model orientation: yaw zero degrees, pitch 
zero degrees; viewed from left side of model; airflow 
left to right 




1 

Fig. 12 - Model orientation: yaw zero degrees, pitch 

zero degrees; viewed from upstream side 
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Fig. 13 - Model orientation: yaw zero degrees, pitch 

4-30 degrees; viewed from left side of model; airflow 
left to right 




Fig. 14 - Model orientation: yaw zero degrees, pitch 

+30 degrees; viewed from upstream side 



61 




Fig. 15 - Model orientation: yaw zero degrees, pitch 

-30 degrees; viewed from left side of model; airflow 
left to right 




Fig. 16 - Model orientation: yaw zero degrees, pitch 

-30 degrees; viewed from upstream side 
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Fig. 17 - Model orientation: yaw 30 degrees, pitch 

zero degrees; viewed from left side of model; airflow 
left to right 




Fig. 18 - Model orientation: yaw 30 degrees, pitch 

zero degrees; viewed from upstream side 
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: !_ J 

Fig. 19 - Model orientation: yaw 30 degrees, pitch 

+30 degrees; viewed from left side of model; airflow 
left to right 




Fig. 20 - Model orientation: yaw 30 degrees, pitch 

+30 degrees; viewed from upstream side 
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j 

Fig* 21 - Model orientation: yaw 30 degrees; pitch 

-30 degrees; viewed from left side of model; airflow 
left to right 





Fig. 22 - Model orientation: yaw 30 degrees, pitch 

-30 degrees; viewed from upstream side 
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— *1 



i I 

Fig, 23 - Model orientation: yaw 60 degrees, pitch 

zero degrees; viewed from left side of model; airflow 
left to right 




Fig. 24 - Model orientation: yaw 60 degrees, pitch 

zero degrees; viewed from upstream side 
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Fig. 25 - Model orientation: yaw 60 degrees, pitch 

+30 degrees; viewed from left side of model; airflow 
left to right 




Fig. 26 - Model orientation: yaw 60 degrees, pitch 

+30 degrees; viewed from upstream side 
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Fig, 27 - Model orientation: yaw 60 degrees; pitch 
-30 degrees; viewed from left side of model; airflow 
left to right 





Fig, 28 - Model orientation: yaw 60 degrees, pitch 
-30 degrees; viewed from upstream side 
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s " 1 

Fig. 29 - Model orientation: yaw 90 degrees, pitch 
zero degrees; viewed from left side of model; airflow 
left to right 




Fig. 30 - Model orientation; yaw 90 degrees, pitch 
zero degrees; viewed from upstream side 
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Fig* 31 - Model orientation: yaw 90 degrees, pitch 

+30 degrees; viewed from left side of model; airflow 
left to right 




Fig, 32 - Model orientation: yaw 90 degrees, pitch 

+30 degrees; viewed from upstream side 
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Fig 0 33 - Model orientation: yaw 90 degrees, pitch 

-30 degrees; viewed from left side of model; airflow 
left to right 




Fig* 34 - Model orientation: yaw 90 degrees, pitch 

-30 degrees; viewed from upstream side 
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Fig. 35 - Model orientation: yaw zero degrees, pitch 

+60 degrees; viewed from left side of model; airflow 
left to right 




Fig. 36 - Model orientation: yaw zero degrees, pitch 

+60 degrees; viewed from upstream side 
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Fig. 37 - Model orientation: yaw zero degrees, pitch 

+90 degrees; viewed from left side of model; airflow 
left to right 




Fig. 38 - Model orientation: yaw zero degrees, pitch 

+90 degrees; viewed from upstream side 
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Fig. 39 - Model orientation: yaw zero degrees, pitch 

-60 degrees; viewed from left side of model; airflow 
left to right 




Fig. 40 - Model orientation: yaw zero degrees, pitch 

-60 degrees; viewed from upstream side 
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p 

Fig. 41 - Model orientation: yaw zero degrees, pitch 

-90 degrees; viewed from left side of model; airflow 
left to right 




1 



Fig. 42 - Model orientation: yaw zero degrees, pitch 

-90 degrees; viewed from upstream side 
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Fig. 43 - Support arrangement used for yaw variation 
in five degree increments; model at zero degrees 
pitch; viewed from side; flow left to right 




Fig. 44 - Same as 43, but view from upstream side 
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FIG. 45 DRAG SUMMARY CONTOUR PLOT 
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FIG. 51 
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FIG. 52 CONVERSION SCALE 
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Fig. 56 - Photograph of the West Coast Research 48 in. 
x 60 in. wind tunnel, showing support beams for wire 
attachment exterior to the test section 




Fig. 57 - The dynamic model in the test section, 
side view 
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Fig. 58 - Photograph of the dynamic model in the 



test section, view from downstream 




Fig, 59 - Close-up of the dynamic model 
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Fig. 60 - Photograph of dynamic model instrumentation 
hardware and bearings 




Fig. 61 - Photograph of the modified potentiometer 
used in dynamic model instrumentation 
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3.0 4.0 5.0 6.0 70 8.0 

DYNAMIC PRESSURE , q ( psf ) 
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DYNAMIC PRESSURE, q (psf) 



FREQUENCY SQUARED, CJ 2 (cps)^ 



1.4 



1.2 



1.0 



0.8 



0.6 



0.4 



0.2 



0 



0 



FIG. 65 




□ MODEL D 

A MODEL E 



2.0 4.0 6.0 8.0 

DYNAMIC PRESSURE, q (psf) 
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FIG. 67 TIME HISTORY OF FREE 
VIBRATION IN YAW 
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